Transient State Work Fluctuation Theorem for a Driven Classical Dissipative System 
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We derive the nonequilibrium transient state work fluctuation theorem and also the Jarzynski 
equality for a classical harmonic oscillator linearly coupled to a harmonic heat bath, which is dragged 
by an external agent. Coupling with the bath makes the dynamics not only dissipative but also non- 
Markovian in general. Since we do not assume anything about the spectral nature of the harmonic 
bath the derivation is not only restricted to the Markovian bath rather it is more general, for a 
non-Markovian bath. 



I. INTRODUCTION 

Although the field of equilibrium thermodynamics and 
equilibrium statistical mechanics are well explored, there 
existed almost no theory for systems arbitrarily far from 
equilibrium until the advent of fluctuation theorems 
(FTs)±'2 i M££ 1 Z i n m id go's. In general, these fluctua- 
tion theorems have provided a general prescription on 
energy exchanges that take place between a system and 
its surroundings under general nonequilibrium conditions 
and explain how macroscopic irreversibility appears nat- 
urally in systems that obey time reversible microscopic 
dynamics. Fluctuation theorems have been cast for var- 
ious nonequilibrium quantities like heat, work, entropy 
production etc and for systems obeying Hamiltonian 8 
as well as stochastic dynamics^ 1 ^. Quantum versions 
of FTs are also know n 11 ' 12 . 

Apart from fluctuation theorems, Jarzynsk i 13 ' 14 had 
also provided a remarkable relation between the work 
done on a system with the equilibrium free energy dif- 
ference. To illustrate the relation, consider a classical 
system in contact with a classical heat bath at temper- 
ature T. Initially the system was in equilibrium with 
the bath and then driven by some external agent (gen- 
eralized force) /. Now the free energy F for the system 
can be calculated by computing the partition function Zf 
when the generalized force is fixed at the value /, since 
Ff = -~{3~ 1 lnZf, where f3 = jr-^, kg is the Boltzman 
constant. Let the system starts in equilibrium at time 
t = specified by / = A and then driven off to a later 
time t = t. If this process is done quasistatically then 
the system remains in equilibrium at each stages of the 
process and also at t = r specified by / = B. Then the 
work done W on the system equals the free energy differ- 
ence, AF = Fb — Fa- On the other hand if the process 
is carried out with a finite rate (i.e. in nonequilibrium 
conditions), W will on average exceed AF. 



(W) > AF 



(1) 



The external agent / is always varied in precisely 
the same manner from A to B. After each realiza- 
tion the work W performed on the system is calcu- 
lated and the distribution function for W, P(W) is con- 
structed. Jarzynski derived the following mathematical 



equality popularly known as the Jarzynski equality (JE) 
or Jarzynski relation where the angular bracket indicates 
the average taken over the distribution function P(W). 



(exp(-pW)} = exp(-/?AF) 



(2) 



Essentially there are two classes of fluctuation theo- 
rems, steady stated and transient fluctuation theorems 
(TFT) 6 . Since the paper deals with the transient state 
fluctuation theorem for work here we do not discuss any- 
thing further on steady state fluctuation theorem. The 
transient state work fluctuation theorem gives the ratio 
of probabilities for the production of positive work to the 
production of negative work as follows, 



P(+W) 
P{-W) 



exp(0W) 



(3) 



Where, W is the work done on the system by an ex- 
ternal agent for an arbitrary time period r. The theorem 
holds for any value of t, provided one starts with the 
system in equilibrium. 

Crooks 1 ^ provided another relation for the dissipative 
work, W dlss , defined as W diss = W — AF. Where AF 
is the free energy difference between the final and the 
initial equilibrium state. The relation is very similar to 
the transient work fluctuation theorem and known as the 
Crooks fluctuation theorem (CFT). 



P F {+W dlss ) 
P R (-W diss ) 



exp([3Wdiss) 



(4) 



Here Pr(— Wdiss) is the probability distribution of neg- 
ative dissipative work done in a time-reversed process. 
Clearly AF = means W = Wdiss- In such a situation 
the transient state fluctuation theorem for work and the 
Crooks fluctuation theorem are equivalent. 

In this paper we derive the transient state fluctuation 
theorem for work for a classical harmonic oscillator cou- 
pled linearly to a harmonic bath. Because of the cou- 
pling to the bath, the system becomes dissipative. We 
start from a Hamiltonian description for the system plus 



2 



the harmonic heat bath and then the system is driven 
by an external agent for a time period of r for a series of 
measurements. We analytically calculate the distribution 
function for work and show that it obeys the transient 
state fluctuation theorem. For our particular choice of 
the external agent the free energy change for the process 
is zero and hence the transient fluctuation theorem and 
the Crooks fluctuation theorem are same. 

In recent past, experiment o 16 ' 17 were carried out to 
verify the FTS. In the experiment by Wang et al. a col- 
loidal particle was trapped using Laser and then dragged 
through a solvent and subsequently the fluctuation the- 
orem was verified. Our Hamiltonian efficiently mod- 
els such situation since the harmonic oscillator in our 
model could be viewed as the colloidal particle in the 
harmonic trap caused by Lasers and the harmonic bath 
as the solvent through which the colloidal particle is 
dragged. People have already verified FTs ^ 18 ' 19 in the 
context of the above mentioned experiment. This mod- 
eling was based on a Langevin description which was 
Markovian. Only very recently Mai ei.ai^S., Speck et.al^l 
and Ohkuma 22 et.al have reported derivations based on 
generalized Langevin equation taking care of the Non- 
Markovian nature of the dynamics. But here we do not 
start from a Langevin description rather we start with 
a Hamiltonian description. Since we couple our system 
Hamiltonian with a set of bath oscillators effectively it 
produces a noise acting on the system which depends on 
bath variable and thus the dynamics becomes stochas- 
tic. Also our derivation does not assume anything about 
the spectral nature of the bath and the particle coupling. 
Hence the results are of very general validity. The pa- 
per is arranged as follows. In the next section we define 
work done on the system. In section III we introduce our 
model. Section IV contains the detailed derivation of the 
fluctuation theorem and in the last section we conclude 
our results. 



II. DEFINITION OF HEAT 



harmonic heat bath, a model for the surrounding solvent. 
In the experiment by Wang, et al the colloidal particle 
was dragged through the solvent by using a Laser trap. 
The particle experiences a harmonic trap whose minima 
moves in time. To model such a situation we assume that 
the particle x is in a harmonic well whose minima is time 
dependent. We thus introduce the Hamiltonia n 23 ! 24 ' 25 

H = H S + H B + h mt (6) 

where 

H s = ^+ k -{x-a{t)f 
2m 2 

and 

^ + ^.=e(!4 (•-$*)") 

Hs is the Hamiltonian for the system, Hb is that for 
the harmonic bath and h int represents the coupling of the 
system with the bath. Here p x and pi are the momenta 
for the particle and the i-th bath coordinate, k is the 
force constant of the optical trap and a(t) is the time 
dependent mean position of the harmonic trap. 

IV. DERIVATION OF TFT FOR WORK AND JE 

The time evolution of the system plus bath is governed 
by the Hamiltonian H. The equations of motion for the 
system and the bath oscillators are 



Px ^ --q^ = -k(x ~ a{t)) -^2~^x + ^2^ (7) 



According to Jarzynskiii, the work done on the system qj_^ 
(described by the Hamiltonian Hs) by an external agent pi = — — — = — Wjft + CiX (8) 

/ acting on the system from t = to t = r is defined as * 

dH p x 

W = / f^gfdt (5) d Px m 



dH 

We use this definition. fa = = p i (10) 

dpi 

III OUR MODEL ^ n our case ^ e ex t erna l agent / is a(t) and hence the 

work done on the system is 

We consider a classical particle of mass m described 
by a positional coordinate x which is linearly coupled to J dH$ f 

a set of harmonic oscillators, each of unit mass described W— a(t)— — dt = / —ka(t)(x(t) — a(t))dt. (11) 

by a positional coordinate qt(i — 1,2, ...A) forming a q ' q 
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Thus within the harmonic Hamiltonian description the 
work done on the system W is linear in x. As our 
Hamiltonian is quadratic and we shall assume equilib- 
rium distributions for the initial distribution for the 
initial conditions, W would have a Gaussian probabil- 
ity distribution^ with a mean (W) and the variance 



*w = (w 2 ) - (wy 

W is 



So the distribution function for 



ka(s) + mx(0) + E c ;<7i( s ) + msx(0) 
*(«) = - 7 1=1 v x 



N c2 

ms 2 + k + J2 



and 



P(W) 



,-(W-(W)) 2 /2a 2 w 



(12) 



W 



with 



(W) = J -ka(t)((x(t)) -a{t))dt (13) 
o 



and 



w 



dh / dt 2 a(t 1 )a(t 2 )C(t 1 ,t 2 ) (14) 



(x(h)) and C(t u t 2 ) 



with Ax(ti) 
{Ax(h)Ax(t 2 )) . 

With the above Gaussian distribution function for 

2W(W) 

P(W) — 2 — 
work it is easy to show that priyb — e " w ■ So in 

order to satisfy the TFT for work it is enough to show 

that a 2 w = *p.. 

Now it is obvious that to calculate the work done, its 

mean and the variance one has to know i as a function 

of time t. To find function of time t we proceed as 

follows. We take a Laplace transform of Eq. ([7]) and Eq. 

©to get. 



J 



{ Pi {0) + sq.jO) + c t x(s)) 



(16) 



Now we substitute g,(s) in Eq. (jTSj) from Eq. (Tl7j)) and 
get x(s) in terms of the initial momenta and position of 
the bath coordinates in time and that of the system itself. 



x(s) = k&{s) + (p a (0) + mx(0)s + g(s))b(s) (17) 
Where 



600 



A" 



JV 



k + ms 2 + E ^ " E J^t^ 

i—1 z 2—1 % 



(18) 



and 



JV 



/(•) = £■ 



(Pi(0) + agi(0)) 
(s 2 +^ 2 ) 



(19) 



Taking the inverse Laplace transform of Eq. (|17[) one 
jets 



at(i) = mx(0)y(i) + mu(0)6(t) + / dt'b(t - t') [ka{t') + £(*')) 

o 



(20) 



with 



conditions. This readily gives 



/v 



£(*) = <?(*) - Y, 4 cos(wit)a;(0) 



and 



y(t) = / dt 7 | -i') + -y4 cos(cjit')b(t - t') ) 

Now Eq. ([2U)l should be consistent with the initial 



my(0) = 1 
m6(0) = 
mb(0) = 1 
my(0) = 



(21) 
(22) 
(23) 
(24) 
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One can substitute x(t) from Eq. (|20|) to Eq. (fTTj) obtained from Eq. (|TT|) by replacing &(f) with its thermal 
to calculate W. Next task is to calculate (W) which is average 



(x(t)) = m (x(0)) y(t) + m («(0)) b(t) + J dt'b{t - t') (fea(t') + (£(£')» 

o 



(25) 



r 



Here the angular bracket indicates a thermal average 
taken in the initial state of the harmonic bath (at t = 0) 
with the shifted canonical equilibrium distribution (since 
we start from an initial equilibrium distribution) given 



J 



by p ~ e-^ (0) . Thus (p x (0)) = 0, (x(0)) = a(0). Now 
to calculate (<?(£)) we proceed as follows. First we take 
the inverse Laplace transform of Eq. (fH?)) to get 



.9(0 = E ^ IpKO) 5 ^ + gi (0)cosM)} 

i=l L J 



(26) 



Next we take the thermal average with respect to the 
initial distribution p ~ e^^ Hi ' to get 



JV c 2 

= E ^x(O) cos(wii) 

i=l ! 



(27) 



as (pi(0)) = 0, ( ft (0)) 



r x(0). Thus (f(f)) 



N 



(9(t)) ~ E cos(ujit)x(0) = 0. The quantity is a 

i=l 1 

Gaussian random force from the bath with the statistical 
properties, (£(*)) = 0, (£(*)£(*')) = /i" 1 ]^ - *'), where 

JV 2 

r(«) = 

Finally we get 



= ma(0)y(t) + / dt' ka(t')b(t - t') (28) 


Next we derive a set of equations in the Laplace and 
the time domain. These are used to get the TFT for 
work. 

In a few steps one can show, msy(s) = 1 — kb(s), 



where y(s) = sb(s) + ^-b(s), T{s) = ~£ ^/ 2 c ;^ r 

i=l * z 

verse Laplace transform of msy(s) = 1 — kb(s) gives 

*(*) = Also si ( s ) = y( s ) - iP( s ) whose 

inverse Laplace gives b(t) = y(t) - i / dt'T(t')b{t - t'). 

o 

So we have the following set of important relations in the 
Laplace and the time domain. 



JV 



In- 



msy(s) = 1 — kb(s) 
sb(s) = y(s) - ±T(s)b(s) 

m = (&) K*) 

b{t) = y{t)-±}dt'T{t>)b{t-t>) 



(29) 



Using the above set of equations one can show that the 
average work done on the system (which is given by Eq. 

(QD) is 



(W) = (y) J dt J dt'ka(t)y(t-t')ka(t') (30) 
o o 

To derive it we proceed as follows. First we replace 
(x(t)) in Eq. JTSJ) from Eq.® to get 



{W) = | _ °m 

r r t 

J dtkrna(0)a(t)y(t) - J dt J dt'ka(t)b(t - t')ka(t') 





and then using Eq. (|29|) followed by an integration by 
parts and using Eq. (f2Tj) one ultimately gets 



(W) = dt dt'ka{t)y{t - t')ka{t') 



o o 
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which is Eq. (|30|). 

In order to evaluate the variance we first have to calcu- 



late the correlation function, C(t\,t 2 ) = (Ax(ti)Ax(t 2 )). 
Now using Eq. ([2T)|) one can show 



J 



C(h,t 2 ) = m 2 {kpy 1 y(ti)y(t 2 ) + m 2 (to/?)" 1 6(ti)6(t 2 ) + /T 1 J dt[ J dt' 2 b(h - *i)6(t 2 - ^TO - f a ) (31) 

o o 

I — 



Where we have used = & y V{t — t'), The above expression for C(t 1} t 2 ) is then put back in 

(x(0) - a(0)) 2 ) = {k(3)' 1 , (v 2 (0)) = {rn(3)-\ E <1- to get 

I 



(w) (l dtka(t)y(t)) + (f ) (} dtka(t)b(t) 



jjjdhf dt 2 J dt[ J dt 2 T(t[ - t! 2 )ka{tx)b(t x - *i)6(t 2 - t' 2 )ka{t 2 



(32) 



r 



Let us define 



ti 



J(h,t 2 ) = J dt\ J dt' 2 r(t; - t' 2 )b(ti - *i)6(ta - Q 



With the help of Laplace and Fourier transform one can 
show (see Appendix for the evaluation of the integral) 



(vn \ ( in \ 

j) y(h-h)- -j tf(*i)i/(* 2 )-m6(ti)6(t 2 ) 

(33) 

When the above expression for J(ti,t 2 ) is plugged into 
Eq. ([311) one § ets 

T T 

a W = (^) y ^ / dt2k ®( tl )y( tl _ t2)ka{t 2 ) 
o o 

In Eq. (|30p . t and i' being dummy variables one can 
change these into t\ and t 2 respectively and rewrite Eq. 



T tl 



(W) = (^) / dti / dt 2 ka(h)y(h-t 2 )ka(t 2 ) (34) 



o o 

Now we interchange t\ and t 2 . 



(W) = (y) / rft 2 / dt x ka{t 2 )y{t 2 -t x )ka{tx) (35) 



Adding Eq. flM)). and Eq. ([35]). Since y(i) is an even 
function of t and also r > t\ , r > t 2 one gets 



2(W) = (^) J dh j dt 2 ka(t 1 )y{t 1 -t 2 )ka{t 2 )=/3a 2 v 
o o 

This shows that the TFT for work is satisfied. Next 
we come to Crooks fluctuation theorem. As we have 
mentioned earlier that TFT for work and CFT becomes 
equivalent when the free energy change for the concerned 
process becomes zero (AF = 0, W = Wdiss) which 
happens in our case. This is because we assumed the 
external agent is the time dependent mean position of 
the harmonic trap and subsequently when one evaluates 
the partition function by performing a Gaussian integral, 
a(t) does not appear in the partition function and thus 
the free energy becomes independent of a(t) so the free 
energy change is zero. Hence the TFT for work reads as 



P(+W) 



= exp(f3W) 



P{-W) 

JE is obtained easily by integrating above equation 

oo 

(exp(-/3W)) = J dWP(+W) exp(-/3VF) 

— OO 

OO 

= J dWP(-W) = 1 



(36) 



o o 
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CONCLUSIONS 
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bath. 
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VII. APPENDIX 

The integral J(ti,t 2 ) is evaluated as follows. First let 
us consider the integral 



J 



/ dt 2 e~ S2t2 J dt' 2 b{t 2 - t' 2 )e- iut 2 
oo 



= b(s 2 ) J dt 2 e' S2t2 e~ luJt2 = b(s 2 )'- 
o 



-(°2+*")*2 , 00 _ b(s 2 ) 

(s 2 +iu)) 10 (s 2 +iLu) 



r 



(37) 



Similarly 



Then with the help of Eq. |j37j and Eq. ([38]) double 
Laplace transform of J(t\,t 2 ) can be written as 



J dt ie ~ Sltl J dt'^ih - t[)e 



(Sl — ILU) 



(38) 



J(si,s 2 )= / dt ie ~ Sltl / dt 2 e~ S2t2 / dt' x / dt^T^ - t' 2 )b{t x - t[)b(t 2 - 1' 2 ) 



o o 



(39) 



We also define 

oo 

r(*i-*2) = -m= I durMeWi-tti and 



Then after some algebraic manipulations Eq. (|39|) be- 
comes 



t \ 1 2 i v 

J(si,s 2 ) = J dt x e~ s ^ J dt 2 e- S2t2 J dt\ J dt' 2 b(t x - t[)b(t 2 - t' 2 ) ( ^- / dujf (w)e^-*i) 

oo 



\T2ir J (si 



b{s 1 )b(s 2 )r(u J ) _ b( S1 )b(s 2 ) 
j)(s 2 +iuj) ( Sl +s 2 )V2 



= f dLJ 

IT J \ S 1- 
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(40) 



- b( Sl )b(s 2 ) f ( T(u) , r( w ) 
J{si,s 2 ) = — I dw 



[si + s 2 )V2tt 



si — iu> s 2 + iu> 



Next consider the integral — i== f du) F ^ which can 

b V27T J Sl— IU 

— OO 

be evaluated as follows 



(41) 
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2tt J Si— ilo v2tt v2tt J C 5 l — ^ v / J V / I — 2-ki J (uj-\-isi) 



Similarly 



(42) 



J dtT(t)Q(t)e- Slt = f dtr{t)e- Slt = f(si) 

-oo 



= / cLj — = f(s 2 ) 
2tt J s 2 + iuj 



Now we use Eq. ((42j) and Eq. (43]) to get 



(43) 



J(si, s 2 ) = y(ti-* a )- f — J y(h)y(t 2 )-mb(h)b{t 2 ) 

(45) 

Taking the double Laplace transformation of the above 
Eq. gives 



^(si,S2) 



b{s 1 ) + b{s 2 ) 
(si + s 2 ) 



(f( Sl ) + f( S2 )) (44) J(ti)fe) = ^ y{tl - t2 )- (Tp) y^MhymbfrMt,) 



The above Eq. is then simplified with the help of the 
Eq. (HH) to get 



which is Eq. 
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